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EXECUTIVE  SUMMARY 


This  work  provides  numerical  modeling  of  the  near  field  of  a  wire  antenna.  The 
conventional  near  field  definition  is  derived  along  with  the  analytic  expression  for  the 
power  density  in  the  near  field.  The  study  is  motivated  by  the  fact  that  the  electric  and 
magnetic  fields  in  the  near  field  region  may  pose  a  radiation  hazard  due  to  higher  than 
expected  field  values.  With  accurate  modeling  of  the  near  fields,  it  will  be  possible  to 
identify  and  predict  such  hazards.  Numerical  results  are  provided  for  the  electric  and 
magnetic  fields  in  the  near  field  region  of  a  thin  wire  antenna.  The  results  obtained  from 
the  model  are  shown  to  compare  very  well  with  data  available  in  the  literature  for  the  thin 
wire  case.  The  variation  of  the  wave  impedance  and  power  density  are  also  studied.  It  is 
determined  that  the  wave  impedance  goes  through  very  sharp  changes  in  the  near  field 
region  and  that  it  approaches  the  free  space  value  as  the  far  field  region  is  approached.  A 
similar  behavior  is  observed  for  the  power  density  as  well.  The  next  step  in  the  study  will 
be  to  consider  the  impact  of  the  near  field  excitation  on  a  dielectric  body.  The  heating 
effect  in  the  dielectric  body  due  to  the  near  field  will  provide  insight  into  the  impact  of 
the  near  field  in  comparison  to  a  plane  wave  excitation.  The  work  can  also  be  extended  to 
consider  the  near  field  of  other  antenna  geometries  of  interest.  For  convenience,  the 
computer  code  developed  in  MATLAB  is  provided  in  the  appendix. 
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1.0  INTRODUCTION 


The  objective  of  this  work  is  to  develop  a  computational  model  to  accurately 
predict  the  near  field  of  an  antenna.  As  a  starting  point  and  for  the  sake  of  simplicity,  a 
thin  wire  antenna  is  considered.  After  modeling  the  near  field  for  this  antenna  and 
understanding  the  mechanics  of  the  associated  near  field,  more  complex  antenna 
configurations,  such  as  a  parabolic  dish  antenna  will  be  attempted  as  a  continuation  of 
this  work.  To  give  a  simple  definition  for  the  near  field,  as  the  name  suggests,  it  is  the 
field  very  close  to  the  antenna.  A  more  rigorous  definition  is  provided  in  the  next  section. 
The  nature  of  the  near  field  is  quite  complex.  As  a  result,  very  few  antenna  geometries 
have  been  studied  in  detail.  Most  often  the  models  make  assumptions  to  simplify  the 
mathematical  analysis  and  modeling.  Thereby,  leaving  the  true  nature  of  the  near  field 
relatively  unknown.  An  accurate  description  of  the  near  field  is  necessary  because  of  its’ 
potential  health  hazard.  This  is  due  to  the  fact  that  in  the  near  field  region,  the  intensity  of 
the  field  and  consequently  the  resulting  power  density  can  exceed  the  recommended 
safety  level.  The  radiation  safety  level  given  by  the  Federal  Communications 
Commission  (FCC)  is  5  mW  I  cm1  for  controlled  exposure  and  lmW/  cm1  for  general 
population  in  the  frequency  range  from  1.5  GHz  to  100  GHz  [1].  The  safety  regulation  by 
OSHA  limits  the  power  density  to  10  mW /cm2  when  spatially  averaged  over  a  6  minute 
period.  Experimental  work  has  demonstrated  that  the  effect  of  radiofrequency  (RF) 
radiation  on  living  tissue  is  primarily  in  producing  heat.  This  heating  is  a  function  of  the 
strength  of  the  RF  radiation  as  well  as  the  duration  of  exposure.  The  depth  of  penetration 
of  the  heat  depends  on  the  frequency  of  the  radiation.  Radiation  at  lower  frequencies  (200 
MHz  to  900  MHz)  penetrates  deeply  whereas  radiation  at  much  higher  frequencies  (1.5 
GHz  and  up)  typically  used  by  radars  produce  heating  effects  at  or  near  the  surface.  The 
human  body  is  capable  of  diffusing  a  portion  of  the  heat  as  circulating  blood  in  the  body 
acts  as  a  coolant.  But,  body  organs  which  have  little  blood  flow,  such  as  eyes  are 
particularly  at  risk  [2].  Other  risks  from  high  levels  of  RF  radiation  include  ignition  of 

1 


Distribution  A:  For  Public  Release 


fuel.  Electromagnetic  energy  is  capable  of  inducing  currents  into  any  metallic  object. 
This  means  that  many  parts  of  an  aircraft  and  refueling  vehicles  may  act  as  receiving 
antennas  and  create  large  enough  power  density  to  cause  a  spark  and  ignite  fuel.  To 
prevent  such  a  potential  hazard,  a  safety  criterion  of  5  W/ cm 2  has  been  established  for 
refueling  operations  [3]. 

2.0  DEFINITION  OF  NEAR  FIEED 

In  this  section,  a  precise  definition  of  the  near  field  in  terms  of  the  distance  from 
the  antenna  is  provided.  Traditionally,  textbooks  in  antenna  theory  define  the  near  field  to 
be  the  region  surrounding  the  antenna  such  that  kr  «  1,  where  k  is  the  propagation 
constant  (or  wave  number  of  the  medium)  and  r  is  the  radial  distance  from  the  origin.  An 
alternate  definition,  which  defines  the  beginning  of  the  far  field  region  (also  known  as 
Fraunhofer  region)  is  derived  as  follows  [4],  Consider  a  filament  of  current  along  the  z- 
axis  and  located  near  the  origin  as  shown  in  the  Fig.  1. 

T-t 

► 

7 

Figure  1:  Field  calculation  for  a  filament  of  current  source  along  the  z-axis. 

The  magnetic  vector  potential,  A , ,  due  to  the  current  distribution  on  the  filament, 
7(z )  ,  is  given  by 

-  jkr 

A.  =  fl(z')— — dz'  (1) 

z  J  4  nR 
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where  R  =  \r  -r  1  represents  the  distance  between  the  observation  or  field  point  ( located 

by  the  position  vector  r  )  and  the  source  point  (located  by  the  position  vector  r').  Since 
A,  is  independent  of  the  angular  variation,  (f) ,  the  vector  potential  can  be  evaluated  in  the 

y-z plane  ,  that  is,  for  (j>  =  90° ,  yielding 

r  =  yy  +  zz,  r'  =  z'z  and  R  =  r  -  f  |  =  *Jy2  +(z-  z')2  .  Converting  rectangular  to  spherical 
coordinates,  y  =  r  sin#,  and  z  =  r  cos  9,  R  can  be  written  as 
R  =  -Jr2  +  z'2  -  2rz'  cos  6* 


1 


R  =  r{  1  +  — (z'2  - 2rz'cos#)} 


1 1/2 


The  above  expression  for  R  can  be  expanded  using  Binomial  theorem  expansion  as 


R  =  r(l  +  (z'2  -  2 rz  cos  0)  -  cos 1  6  +  terms  of  order  f  \  ]  + . ..) 


=  r-z'cos6*  +  — z'2(l-cos2  0)  + ... 

2  r 

=  r-z'cos6l  +  —  z'2  sin2  9  +  ...  (3) 

2  r 

In  the  integral  given  by  Equation  (1),  the  factor,  R,  in  the  denominator  can  be 
approximated  by  R  «  r .  However,  in  the  phase  term,  more  accuracy  is  required  and  the 
proper  approximation  is 

R»r-z'  cos  0  (4) 

This  approximation  is  illustrated  in  Figure  2  as  rays  are  drawn  from  each  point  on  the 
source  as  parallel  lines.  The  parallel  ray  assumption  is  exact  only  when  the  observation 
point  is  at  infinity,  but  it  is  a  good  approximation  in  the  far  field.  The  definition  of  the 
distance  from  the  source  where  the  far  field  begins  is  taken  to  be  where  the  parallel  ray 
approximation  begins  to  breakdown.  The  distance  where  the  far  field  begins,  rff  is  the 

z' 

value  of  r  for  which  the  path  length  deviation  due  to  neglecting  the  term  — sin2  0  in 

2  r 
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Equation  (3)  is  —  .  This  corresponds  to  a  phase  error  of^-—  =  —  radians  .  If  2)  is  the 
16  F  F  /l  16  8 


length  of  the  line  source,  rff  is  found  by  equating  the  maximum  value  of  —sin"  6  equal 


2  r 


to  — ,  that  is, 
16 


<2 


2  r 


-sin2  0 


ff 


A 

16 


This  occurs  for  z'  =  D / 2  and  6  =  90° ,  yielding 
2  D2 


rff  = 


A 


(5) 


With  the  far  field  distance  identified,  the  distance  to  the  near  field  is  taken  to  be  half  of 
this  distance.  Therefore,  the  extent  of  the  near  field,  rnf ,  is  given  by 


V 


D 2 

"T 


(6) 


The  region  between 


D2 

T 


4  2Z)2 
and  - 

A 


is  termed  as  Fresnel  region  or  intermediate  near  field 


region,  as  illustrated  in  Figure  3. 


Figure  2:  Parallel  ray  approximation  for  far  field  calculation  of  a  filament  of  current. 
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Figure  3:  Near  field  and  far  field  regions. 

3.0  UPPER  BOUND  ON  THE  POWER  DENSITY  IN  THE  NEAR  FIELD 

As  the  power  density  in  the  near  field  defined  by  W ,  can  exceed  the  power  density 
at  the  antenna  aperture  defined  byW0,  it  is  instructive  to  know  the  upper  bound  on  W 

that  has  been  established  in  the  literature  [2].  Consider  an  electromagnetic  (EM)  plane 
wave  impinging  on  a  totally  absorbing  body  of  area  A  .  If  the  power  density  of  the  EM 
wave  is  W ,  then  the  power  absorbed  P  is  given  by 
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P  =  WA  (7) 

For  an  isotropic  radiator  in  free  space,  radiating  a  total  average  power  P  in  all  directions 
equally,  then  the  power  density  on  the  surface  of  a  concentric  sphere  of  radius  r  is  given 
by 


P  P 
W  =  -  = 


(8) 


A  4  nr2 

In  case  the  radiator  is  not  isotropic  because  it  radiates  with  a  directive  gain  G{6,  (f>)  in  a 
given  direction,  then  the  power  density  in  the  far  field  at  a  distance  r  will  be 
GP 


W  = 


4n  r 


(9) 


A  100%  ground  reflection  doubles  the  strength  of  the  electric  field  and  hence  the  power 
density  gets  quadrupled,  That  is 
GP 


W  = 


n  r 


The  antenna  gain  G  can  be  written  in  terms  of  the  antenna  area  as 
4  n  A 


G  = 


/l2 


(10) 


(ID 


Substituting  G  from  Eq  (1 1)  into  Eq  (9)  gives  the  far  field  power  density 
W=  AP  _W0A2 


A  V  A2r2 

where  W0  =  PI  A  is  the  power  density  at  the  aperture.  The  above  equation  can  be  written 
as 


W  A2 


2  2 


W0  Azr 


(12) 


Similarly,  substituting  G  from  Eq  (11)  into  Eq  (10)  gives  the  power  density  for 
the  100%  reflection  from  the  ground: 
w_  4AP  _  4  W0A2 


i2  2  i2  2 

A  r  A  r 


or 
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The  expression  in  Eq  (13)  is  valid  in  the  far  field  region,  but  for  a  uniformly  illuminated 


round  aperture  and  in  the  limit 


field  region  providing 


f  A 2  ^ 

v^Vy 


«  1  ,  the  expression  can  be  applicable  in  the  near 


W  .  .  2 
—  =  4  sin 
Wn 


A  r 


(14) 


o 

which  has  a  maximum  value  of  4,  thereby,  providing  the  approximation  for  the  near  field 
as 


W  =  4W0  (15) 

The  power  density  of  a  reflector  antenna  can  be  shown  to  have  a  maximum  envelop  along 
the  central  axis  of  the  antenna  which  is  6  dB  (10  log4)  above  the  average  power  density 
on  the  aperture  up  to  a  characteristic  distance  of  D2  /  4 A  and  falls  uniformly  from  that 
point  on  to  the  far  field  [5],  For  instance,  for  a  transmit  power  of  about  3.6  kW  from  a  12 
meter  diameter  aperture,  the  average  power  density  is  about  3.2  mW /cm2,  but  could 
reach  12.8  mW  / cm 2  along  the  central  axis. 


4.0  NEAR  FIELD  OF  A  WIRE  ANTENNA 


This  section  provides  a  mathematical  formulation  for  computing  the  near  field  for 
a  thin  wire  antenna.  A  formulation  for  solving  an  integral  equation  for  the  current 
distribution  on  the  antenna  element  is  presented.  The  resulting  current  distribution  is  then 
utilized  in  computing  the  electric  and  magnetic  fields  in  the  near  field  region.  Consider  an 
arbitrarily  oriented  as  shown  in  Figure  4.  The  wire  is  assumed  to  be  perfectly  conducting. 
For  thin  wires,  it  is  assumed  that  the  wire  radius  is  small  compared  to  the  wavelength  and 
the  length  of  the  wire.  The  thin  wire  assumption  allows  the  consideration  of  only  axially 
directed  currents.  Therefore,  there  is  no  azimuth  or  circumferentially  directed  component 
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of  the  current.  The  magnetic  vector  potential  A  and  the  electric  scalar  potential  ®  are 
given  by 


A  =  — -  [  I  (s')  s  (5')  k(s  -  s’)  ds' 

c 

and 


(16) 


®  = 


- —  [  q(s')  k(s  -  s')  ds' 
4/T£J 

c 

respectively,  where 

1  r*  e~JkR 

k(s-s')  =  —  - ds’ 

2 R 


(17) 


R  = 


(s-s')2  +4o2sin2  — 
2 


q{s)  =  —^~i(s) 
jco  ds 


(18) 

(19) 

(20) 


In  the  above  equations,  I(s)  is  the  current  along  the  wire  axis  and  q(s)  is  the  associated 
charge. 


Figure  4:  An  arbitrarily  oriented  wire  of  radius  2a. 
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The  wire  is  illuminated  with  an  incident  field  E"'c  and  as  a  result  produces  a  scattered 
field  Es  arising  from  the  induced  current.  Applying  the  boundary  condition  on  the  surface 
of  a  perfect  conductor  that  the  tangential  component  of  the  total  electric  field 
must  vanish  yields: 

(fr"c+ff')-.?  =  0  (21) 

The  scattered  field  can  be  written  in  terms  of  vector  and  scalar  potentials  as 

Es  =  -jcoA-VO  (22) 

Substituting  Eq.  (22)  into  Eq.  (2 1 )  gives 

-Einc -s  =  -jcoA-s-VO-s  (23) 

The  substitution  of  Equations  (16)  and  (17)  into  Eq  (23)  yields  an  integral  equation  which 
can  be  solved  numerically  via  the  method  of  moments  [6].  The  wire  is  subdivided  into 
N  + 1  segments,  as  shown  in  Figure  5. 


Sn-1/2 

n-1 


t 

Q 

Figure  5:  Arbitrarily  oriented  wire  antenna  with  segmentation  scheme. 
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The  current  and  charge  are  expanded  in  terms  of  pulse  basis  functions.  Note  that  the 
current  and  charge  pulses  are  off  by  half  a  pulse  width.  For  the  current  there  is  a  half 
pulse  of  zero  width  on  each  end  of  the  wire.  This  is  not  the  case  for  the  charge  pulses 
which  span  the  entire  wire  length.  The  current  expansion  can  be  written  as 

I(s)  =  f ]l„Pn(s)  (24) 

n= 1 


where  the  pulse  function  is  defined  as 

h 

1 0,  otherwise 


P,Xs)  = 


Sn- 1/2  <  S  <  Sn+ 1/2 


(25) 


The  variable  sn  represents  the  arc  length  to  the  end  of  segment  n  measured  from  one  end 
of  the  wire.  The  position  vector  rn  locates  the  end  of  wire  segment  n  .  The  unit  vector 
sn+ 1/2  is  parallel  to  the  wire  axis  on  segments  +  1  and  defined  as: 


(26) 

K+ 1  ~rn\ 

Testing  the  integral  equation  in  (23)  with  pulse  functions  yields  a  system  of  linear 
equations  defined  by 

[z.T.HrJ  (27) 

where 


Z . = 


-1 


juntos 


k  (w-VwHi, 

\T/  \T/ 

Am-l/2,«,n+l  A  m-l/2,«-l 


VT>  S  VT/ 

/i +1/2  t»!,«,«4-1/2  t  ^;i-1/2  t  nun-112 


VT/  VT/ 

\  m+l/2,n,n+l  .  m+\ /2,n-l,n 

n  / 


s„  -s 


n  — 1 


+ 


s„  -s 


n- 1 


K=Einc(smHrm +V1-7m_vi) 


'F  =  f  k(s  -  s')  ds' 

m,u,v  J  \  m  ' 


rn+l+rn 


'n  4-1/2 
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\Zmn]  is  a  square  matrix  with  m,n  =  1,2 ,...N  .[ln]  and  \Vm ) are  column  vectors  of  length 
N  .  The  column  vector  [Vm  ]  represents  an  applied  voltage  at  the  location  of  the  antenna 
feed  point.  For  an  antenna  that  is  fed  at  the  center  with  V  volts,  all  the  elements  of  [Vm] 

are  set  equal  to  zero  except  the  location  that  matches  with  the  feed  segment  (the  middle 
segment  for  the  center-fed  case). 


Near  Field  Computation 

The  electric  near  field  and  the  magnetic  near  field  can  be  calculated  using  the 
current  distribution  obtained  from  the  numerical  procedure  outlined  earlier.  The  method 
is  based  on  the  work  by  Adams  [7],  [8].  To  compute  the  electric  field  at  a  given  point  in 
the  near  field,  a  small  (testing)  thin  wire  dipole  of  length  A /  is  placed  at  the  point  with  its 
axis  parallel  to  the  vector  component  of  interest.  An  additional  expansion  function  is 
assumed  over  this  test  dipole  so  that  the  total  number  of  unknowns  is  ( N  + 1) .  If  the  test 

dipole  is  open  circuited,  then  IN+1  =  0  and  the  open  circuit  voltage  is  given  by 

<28> 


The  electric  field  along  the  direction  /  is  now  given  by 

E-i  =  — — (29) 
A  / 

This  process  of  placing  the  test  dipole  at  the  location  of  the  near  field  point  is  repeated  for 
each  component  (EX,EV,EZ)  as  well  as  for  each  near  field  point.  The  average  value  of  the 


electric  field  is  given  by 


E..g  = 


t(£,2+£,2+£,2) 


1/2 


(30) 


and  this  is  a  conservative  estimate  of  the  maximum  value.  The  maximum  or  peak  value  is 
given  by 


E peak  ~ 


t(£/+£,2+£;2)  +  t(/l2+£2) 


-il/2 


2  \l/2 


(31) 
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where 


A  =  E2X  cos(2<9v)  +  Ey  cos(2<9r)  +  E]  cos(26*.) 

B  =  E;  sin(2<9v)  +  E;  sin(20v)  +  E:  sin(20z) 

and  0x,0v,0_arQ  the  phase  angles  for  the  x,y,z  components.  The  near  magnetic  field  is 
calculated  in  a  similar  way  as  the  electric  field  in  terms  of  placing  the  test  dipole  and  the 
computation  is  done  using  the  vector  potential.  The  magnetic  field  H  can  be  computed 
from  the  vector  potential  A  as 


H  =  —V  xA  =  — 


B 


B 


dA^_dAy\ 
dry  dz 


1 

X  H - 

B 


(  dA „  dA,^ i 
dx  j 


dz 


y  +  ■ 


B 


dAy  dAx\ 
dx  dy 


(32) 


Where  x,y,z  are  the  unit  vectors  along  the  x,y,z  coordinate  axes,  respectively.  If  the  test 
dipole  is  short  enough  such  that  the  fields  vary  smoothly  over  its’  length  then  the  partial 
derivatives  in  Equation  (32)  can  be  replaced  by  finite  differences: 


H  = 


_1_ 

M 


AA_  A A 


Ay  A z 


x  +  ■ 


_1_ 

B 


A Ax  A A 

V  A z 


Ax 


y  +  ■ 


_i_ 

B 


Ax  Ay 


(33) 


For  example,  A Az  /  Ay  is  the  change  in  the  z  component  of  the  vector  potential  along  a 
y  directed  test  dipole  of  length  Ay  located  at  the  near  field  point.  The  average  and  peak 
values  of  the  magnetic  field  are  computed  using  similar  expressions  as  given  for  the 
electric  field  in  Equations  (30)  and  (3 1)  with  E  replaced  by  //  . 


5.0  RESULTS  AND  DISCUSSION 


In  this  section,  numerical  results  on  the  near  fields  for  a  thin  wire  antenna  are 
presented.  The  wire  antenna  is  symmetrically  located  with  its  orientation  in  the 
z  direction  atx  =  0 ,  y  =  0 ,  as  shown  in  Figure  6,  and  is  center  fed  by  a  voltage  of  1  volt 
at  a  frequency  of  300  MHz.  The  antenna  is  sub-divided  into  N  + 1  segments,  as  shown  in 
Figure  5,  where  represents  the  number  of  pulse  basis  functions  employed  for  the 
current  distribution  on  the  wire.  Figure  7  shows  the  current  distribution  on  a  center-fed, 
half-wavelength  dipole  antenna  of  radius  0.005  A  ,  number  of  pulses  (unknowns)  /V  =  31. 
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The  current  distribution  is  utilized  in  calculating  the  near  fields  using  the  formulation 
outlined  in  the  previous  section.  A  test  dipole  of  length  A /  =  0.00  l/l  in  Equation  (29)  is 
used.  Figures  8  and  9  show  the  computed  near  electric  fields,  Ez  (z-component  of  Electric 
Field)  andE^  ( p  -component  of  Electric  Field),  respectively,  for  the  half- wavelength 
dipole  as  a  function  of  z/X  for  different  values  of  radial  distance  pi X  .  It  can  be  seen 
from  these  figures  that  the  radial  component  E  has  peaks  at  the  gap  ( z  =  0)  and  much 
higher  peaks  at  the  end  of  the  dipole (z/ X  =  0.25) .  The  tangential  component  of  the 
electric  field  E z  is  relatively  small  except  in  the  gap  region  and  near  the  end  of  the  dipole 
antenna.  Figure  10  shows  the  associated  0  component  of  the  magnetic  field  for  this 

antenna  configuration.  Figures  11  and  12  show  the  computed  near  electric 
fields,  if.  and E  ,  respectively,  for  a  0.4484 X  dipole  as  a  function  of  z/X  for  different 

values  of  radial  distance  pi X  .  A  copy  of  the  computer  code  (developed  in  MATFAB) 
used  to  generate  the  result  in  Figure  1 1  is  provided  in  Appendix  I.  Figures  13  and  14,  the 
results  obtained  in  this  work  for  E,  and  E ,  respectively,  for  the  0.4484  2  dipole 

for/?/ X  =0.03  are  compared  with  those  obtained  by  Adams  [8].  The  comparison  with 
published  results  is  good,  thereby,  providing  validity  to  the  near  field  computation.  The 
results  for  the  0.5  X  dipole  presented  in  Figures  8  and  9  also  followed  closely  with 
published  work  of  Adams  [7].  The  impedance  and  power  density  for  the  0.5  X  dipole  are 
shown  in  Figures  15  and  16,  respectively.  The  wave  impedance  in  Figure  15  approaches 
the  free  space  value  of  377  ohms  as  p  approaches  the  far  field  distance,  rff ,  of  0.5  X 

(rfl  =2D2  / X ,  where  D  =  XI 2  is  the  dimension  of  the  antenna).  It  is  interesting  to  notice 

that  the  impedance  graph  in  Figure  15,  for  the  case  of  zl X  =  0.25,  does  not  seem  to 
follow  the  trend  of  lower  impedance  values  with  increasing  values  of  z/X  .  A  reasonable 
explanation  for  this  case  is  that  the  electric  field  close  to  the  tip  of  the  antenna  is  very 
large,  as  seen  from  Figure  8,  giving  rise  to  large  values  of  wave  impedance. 
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6.0  CONCLUSIONS 


The  near  field  and  far  field  regions  are  identified  in  terms  of  the  antenna  dimension.  The 
near  field  of  a  thin  wire  antenna  is  computed  for  a  half  wave  dipole.  The  method  is 
applicable  to  any  thin  wire  antenna  as  long  as  the  wire  radius  is  small  compared  to  the 
wavelength,  typically  for  radius  a  <  0.00  IT  .  It  is  noted  that  the  electric  field  variation  in 
the  near  field  is  quite  different  than  the  traditional  far  field  behavior.  The  electric  field 
amplitude  has  peaks  at  certain  locations  which  give  rise  to  larger  power  density  than 
conventionally  expected.  As  the  electric  field  undergoes  through  sharp  peaks  in  the  near 
field,  the  power  density  also  shows  a  peaking  behavior  before  following  the  far  field 
trend.  The  wave  impedance  in  the  near  field  also  undergoes  sharp  variations  with  distance 
from  the  antenna  and  approaches  the  free  space  value  as  the  far  field  region  is 
approached.  The  numerical  model  accurately  characterizes  the  near  field  as  the  computed 
results  are  compared  with  results  given  in  the  literature.  The  model  will  therefore  provide 
accurate  estimation  of  near  field  so  that  radiation  safety  guidelines  can  be  appropriately 
defined. 

7.0  RECOMMENDATIONS 

The  present  work  modeled  the  near  field  of  a  thin  wire  antenna.  The  near  field  data  from 
this  work  can  now  be  imported  as  an  input  to  the  problem  of  computing  the  heating 
effects  in  a  dielectric  model  of  the  human  body.  It  is  also  recommended  to  model  the  near 
field  of  other  antenna  geometries  such  as  a  parabolic  reflector  antenna. 
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Figure  6:  Wire  antenna  oriented  symmetrically  along  the  z-axis  with  field  point 
(x,y,z). 
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Figure  7:  Current  distribution  on  a  center-fed,  half-wavelength  dipole  (length 
radius  =  0.005/1 ,  frequency  =  300  MHz). 
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Electric  Field  in  Volts/meter 
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z-comp.  of  Electric  Field  for  0.5  lambda  dipole 


12 

10 

8 

6 

4 

2 

0 

0  0.05  0.1  0.15  0.2  0.25  0.3  0.35 

z/lambda 

Figure  8:  Near  field  E.  (z-component  of  Electric  Field)  of  center-fed  dipole  antenna 
(length  =  0.5  A  and  radius  =  0.005  A  ). 
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Electric  Field  in  Vols/meter 
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Figure  9:  Near  field  E  (p  -component  of  Electric  Field)  of  center- fed  dipole  antenna 
(length  =  0.5  A  and  radius  =  0.005  A  ). 
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Magnetic  Field  in  Amps/meter 


Magnetic  Field  (phi-component) 


0.045 

0.04 
0.035 

0.03 

0.025 

0.02 

0.015 

0.01 

0.005 

0 

0  0.05  0.1  0.15  0.2  0.25 

z/lambda 

Figure  10:  Near  field  ( (j)  -component  of  Magnetic  Field)  of  center-fed  dipole  antenna 
(length  =  0.5  X  and  radius  =  0.005  X  ). 
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Electric  Field  in  Volts/meter 


20 


z-comp.  of  Electric  Field  for  0.4484  lambda  dipole 
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Figure  11:  Near  field  E.{  z  -component  of  Electric  Field)  of  center-fed  dipole  antenna 
(length  =  0.4484  A  and  radius  =  0.00503  A  ). 
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Figure  12:  Near  field  E  (p  -component  of  Electric  Field)  of  center-fed  dipole  antenna 
(length  =  0.4484  A  and  radius  =  0.00503  A  ). 
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z  comp  of  Electric  Field  for  0.4484  lambda  dipole 
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Figure  13:  Near  field  z  -component  of  Electric  Field)  of  center-fed  dipole  antenna 
(length  =  0.4484  A  and  radius  =  0.00503  A  ). 
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Figure  14:  Near  field  E  (p  -component  of  Electric  Field)  of  center-fed  dipole  antenna 
(length  =  0.4484  A  and  radius  =  0.00503  A  ). 
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Figure  15:  Wave  Impedance  for  center-fed  dipole  (length  =  0.5  A  and  radius  =  0.005  A  ). 
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Figure  16:  Power  density  for  center-fed  dipole  (length  =  0.5  A  and  radius  =  0.005  A, ). 
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APPENDIX  I :  MATLAB  CODE  FOR  ANTENNA  NEAR  FIELD 


%  z-comp  of  Electric  Field  (Near  Field  Calculation  for  a  Wire 
%Antenna)  for  near  field  calculation  August,  2007 

clear  all 

%********************  input  ******************************* 

%  Variables  with  asterik  (*)  on  the  comment  on  the  right  side  are  INPUT 
%  VARIABLES 


nunkns=31;  %  ********  Number  of  Unknowns  (for  current)  on 

the  Antenna 
nunkns2=nunkns+2 ; 

freq=300*10A6;  %  *******************  Frequency  (Hz) 

vel=3*10A8;  %  velocity  of  light  in  free  space 

omega=2*pi*freq;  %  angular  frequency  (rad/sec) 
wk=omega/vel;  %  propagation  constant  (2*pi/lambda) 

wavelength=vel/freq;  %  wavelength  (m) 

length=0 . 2242*wavelength;  %  ******************  Half  dipole  length 
(in  lambda) 
dl ength=2* length; 

rad=0 . 005*wavelength;  %  ********************  Radius  of  Dipole  (in 

lambda) 
waveimp=377 ; 

(ohms ) 

eps=l/ (36*pi*10A9)  ; 
zg=0  ; 

for  Center-fed) 
delta=2*length/ (nunkns+1) ; 


%  **********  Wave  Impedance  in  free  space 

%  **********  Free  space  permittivity  (F/m) 

%  *********  Feed  point  of  the  antenna  (zg=0 

%  Antenna  segment  length  (2*L/ (N+l)) 


%  Grid  points  (x,y,z)  where  near  field  needs  to  be  computed 
dx=0  ; 
dy=  0  ; 
dz=0 . 001 ; 
ndx=l ; 
ndy=l ; 
ndz=300 ; 

xstart=0 . 0*wavelength; 
ystart_array ( 1 ) =0 . 01 ; 
ystart_array (2) =0 . 03; 
ystart_array ( 3 ) =0 . 05 ; 
zstart=0 . 001*wavelength; 


S-  -k  ~k  -k  ~k  -k  ~k 
o 

o, 

o 


For  Antenna  case,  the 
Forcing  Function 


Dipole  is  Center-fed  with 
"vmvector"  is  Unity  at 


1  Volt  ****** 
feedpoint 


for  p=l:nunkns 
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matchpnt= ( (nunkns-1 )  /  2  )  +1  ; 
if  p==matchpnt 

vmvector (p, 1 ) =1 . 0  ;  %  *******  Antenna  feed  voltage  is  1  volt 

else 

vmvector (p,  1 ) =0 ; 
end 

end 


%  **************************  End  of  Input  Data******************** 

%  The  array  rx,ry,rz  gives  the  x,y,z  components  of  wire  segment  end 
points 

for  nn=l:nunkns2 
rx (nn) =  0.0; 
ry (nn) =  0.0; 

rz (nn)  =-length+delta* (nn-1 ) ; 

end 

%  Calculation  of  the  Impedance  Matrix  "zmatrix" 

factor=-l/ ( j *4*pi*omega*eps) ; 
for  m=l rnunkns 
mp 1 =m+ 1 ; 

[rxl, ryl, rzl] =  rmhvector (rx, ry, rz,mpl) ; 

[ rx2 , ry2 , rz2 ] =  rmhvector ( rx, ry, rz , m) ; 
dif fx=rxl-rx2 ; 
dif fy=ryl-ry2 ; 
dif f z=rzl-rz2 ; 
for  n=l : nunkns 
npl=n+l ; 

%  Contribution  due  to  vector  potential 

psil=vecpot ( rx, ry, rz , wk, rad, mpl ,  npl ,  0.5); 
psi2=vecpot (rx, ry, rz , wk, rad, mpl ,  npl ,  -  0 . 5 )  ; 

%Contribution  due  to  scalar  potential 

psi3=scalarpot ( rx, ry, rz , wk, rad, mpl ,+0.5, npl , npl+1 ) ; 
psi4=scalarpot (rx, ry, rz , wk, rad, mpl ,+0.5,  npl-1 ,  npl )  ; 
psi5=scalarpot ( rx, ry, rz , wk, rad, mpl ,-0.5, npl , npl+1 ) ; 
psi6=scalarpot ( rx, ry, rz , wk, rad, mpl , -0 . 5 , npl-1 , npl ) ; 

%  S  unit  vectors 

[sxl, syl, szl, rmag] =sunit ( rx, ry, rz , npl ) ; 

[sx2,  sy2,  sz2, rmag] =sunit (rx, ry, rz, n) ; 

%  Dot  products 

dotl=psil* (dif fx*sxl+dif fy*syl+dif f z*szl ) ; 


28 


Distribution  A:  For  Public  Release 


dot2=psi2* (dif fx*sx2+dif fy*sy2+dif f z*sz2 ) ; 
dotprod=wkA2* (dotl+dot2) ; 

zmatrix (m, n) =factor* (dotprod-psi3/delta+psi4/delta+psi5/delta- 
psi6/delta) ; 
end 

end 


%  SOLUTION  BY  MATRIX  INVERSION 

solvector=inv ( zmatrix) *vmvector;  %  solvector  is  the  solution  vector 
%with  current  on  the  antenna 


%  Plot  the  Current  Distribution  on  the  Wire  Antenna 

rsolvec=real (solvector) ; 
isolvec=imag ( solvector ) ; 
for  ip=l:nunkns 

rrealpart (ip) =rsolvec (ip) ; 
iimagpart (ip) =isolvec (ip) ; 
zpart (ip) =rz (ip) ; 

end 

for  ipp=l : nunkns2 
if  ipp==l 

realpart (ipp) =0; 
imagpart (ipp) =0; 
elseif  ipp==nunkns2 
realpart (ipp) =0; 
imagpart (ipp) =0; 

else 

realpart (ipp) =rrealpart (ipp-1)  ; 
imagpart (ipp) =iimagpart (ipp-1) ; 

end 

end 


%  If  following  three  lines  are  commented  out  The  Antenna  Current  is 
%not  plotted 

%  plot (rz, realpart, ' x ' , rz, imagpart, 'o') 

%  xlabel (' Length  in  meters') 

%  ylabel ( ' Induced  Current  in  AMPS ' ) 


Near  Field  Computations 
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mpl=nunkns2  +  l  ; 
for  kk=l:3 
nfield  points=0; 
ystart=ystart_array ( kk)  ; 
for  i=l : ndx 

xd (i) =xstart+ ( i — 1 ) *dx; 
xdi=xd ( i ) ; 
for  j  =1 : ndy 

yd ( j ) =ystart+ ( j -1) *dy; 
yd j  =yd ( j ) ; 
for  k=l : ndz 

nfield  points=nf ield  points+1; 
zd (k) =zstart+ (k-1) *dz; 
zdk=zd ( k) ; 

%  ef ieldx, ef ieldy  and  efieldz  are  the  x,y,z  components  of  the 
%Electric  Field 

%  ****  note  *****  in  this  example  only  the  z-comp  of  E-field  is 
%being  computed 

%  therefore,  efieldx  and  efieldy  are  commented  out 


%ef ieldx=ef ieldnear (rx, ry, rz , xdi , ydj , zdk, solvector , wk, rad, factor,  delta, 
%wavelength, 1,0,0,  nunkns )  ; 

%efieldy=ef ieldnear (rx, ry, rz, xdi, ydj , zdk, solvector, wk, rad, factor, delta, 
%wavelength, 0,1,0,  nunkns )  ; 

ef ieldz=ef ieldnear ( rx, ry, rz , xdi , ydj , zdk, solvector,  wk,  rad,  factor,  delta,  w 
avelength, 0,0,1, nunkns )  ; 


%  nearfield ( kk, nfield  points) =abs (efieldx) ; 

%field  is  needed 

%  near f ield ( kk, nf ield_points ) =abs (efieldy) ; 

%field  is  needed 

nearfield (kk, nfield  points) =abs (efieldz) ; 

%field  is  needed 


if 

x-comp 

of 

near 

if 

y-comp 

of 

near 

if 

z-comp 

of 

near 

end  %  loop  over  k  (z  values) 
end  %  loop  for  j  (y  values) 
end  %  loop  over  i  (x  values) 
end  %  loop  over  kk 


%  Plot  of  Electric  Field  vs  z/lambda  for  three  values  of  rho/lambda 


plot ( zd ( : 
hold  on 

) , nearfield ( 1 , : 

)  ,  'b- 

plot ( zd ( : 

)  ,  nearfield (2 ,  : 

)  ,  'r- 

plot ( zd ( : 

) , nearfield ( 3 , : 

) ,  'g- 
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title (' z-comp .  of  Electric  Field  for  0.4484  lambda  dipole') 
xlabel ( ' z/ lambda ' ) 

ylabel (' Electric  Field  in  Volts/meter') 

gtext ( ' rho/ lambda=0 . 01 ' ) 

gtext ('0.03') 

gtext ('0.05') 

hold  off 


9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'S-9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'S-9'9' 

oooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooo 


%%%  FUNCTION  USED  IN  THE  MAIN  PROGRAM 
%%&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& 


&&&&&&&&&&  FUNCTION  "rmhvector 


rt 


&&&&&&&&&&&&&&& 


function [ rmhx, rmhy, rmhz ]  =  rmhvector (rx, ry, rz, p) 

%Gives  the  x,y,z  components  of  vector  locating  the  mid-point  of  a 
%segment 

rmhx= (rx (p+1) +rx (p) ) / 2; 
rmhy= (ry (p+1) +ry (p) ) /2; 
rmhz= ( rz (p+1 ) +rz (p) ) /2 ; 

%  FUNCTION  "sunit" 


function  [ snx, sny, snz , rmag]  =  sunit (rx, ry, rz, p) 

%Gives  x,y,z  components  of  the  unit  vector, s  (p+1/2),  on  segment  p  to 
%p+l 

rxl=rx (p+1 ) ; 
ryl=ry (p+1) ; 
rzl=rz (p+1 ) ; 
rx2=rx (p) ; 
ry2=ry (p) ; 
rz2=rz (p) ; 

rmag=sqrt ( (rxl-rx2) A2+ (ryl-ry2) A2+ (rzl-rz2)  A2)  ; 
snx= (rxl-rx2 ) / rmag; 
sny= (ryl-ry2 ) / rmag; 
snz= (rzl-rz2) / rmag; 

%&&&&&&&&  FUNCTION  “scalarfun”  &&&&&&&&&&&&&&& 


function  psi=scalarfun (rx, ry, rz, wk, rad, m, n, del) 

%psi  function  calculates  the  vector  potential  for  H-field 

rmx=rx (m) ; 
rmy=ry (m) ; 
rmz=rz (m) ; 
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[rmx, rmy, rmz] =rmhvector (rx, ry, rz,m) ; 


if  del==0 . 5 

rnx=rx (n) ; 
rny=ry (n) ; 
rnz=rz (n) ; 

[sxl,syl,szl,  rmag] =sunit (rx, ry, rz, n) ; 

end 

if  del==-0 . 5 

[rnx,  rny, rnz] =rmhvector (rx, ry, rz,  n-1)  ; 

[sxl,syl,szl, rmag] =sunit (rx, ry, rz, n-1) ; 

end 

delta2=rmag/ 2.0; 
if  m==n 

F=@ (s)  exp (-j *wk*sqrt ( (rmx-rnx-s*sxl) . A2+ (rmy-rny-s*syl ) .  A2  + (rmz 
rnz-s*szl) . A2+rad. A2) ) . / . . . 

sqrt ( (rmx-rnx-s*sxl) . A2+ (rmy-rny-s*syl ) . A2+ (rmz-rnz- 
s*szl)  . A2+rad. A2) ; 

psi=quad (F, 0.0, delta2) ; 

%  Self  term  calculation 

else 

%  Non-self  term  calculation  (reduced  kernel) 

F=@ (s)  exp (-j  *wk*sqrt ( (rmx-rnx-s*sxl)  . A2+ (rmy-rny-s*syl )  . A2+ (rmz 
rnz-s*szl) . A2+rad. A2) ) . / . . . 

sqrt ( (rmx-rnx-s*sxl) . A2+ (rmy-rny-s*syl) . A2+ (rmz-rnz- 
s*szl)  . A2+rad. A2) ; 

psi=quad (F, 0.0, delta2) ; 

end 

%  &&&&&&&&&&  function  "scalarpot^  &&&&&&&&&&&&&&&&&&&&&&& 


function  psi=scalarpot (rx, ry, rz, wk, rad, m, del, n, q) 

%psi  function  due  to  scalar  potential 

rnx=rx (n) ; 

rny=ry (n) ; 

rnz=rz (n) ; 

rnx2=rx (q) ; 
rny2=ry (q) ; 
rnz2=rz (q) ; 

[sxl,syl,szl, rmag] =sunit (rx, ry, rz,  n)  ; 
if  del==0 . 5 

[ rmx, rmy, rmz ] =rmhvector ( rx, ry, rz , m) ; 

end 

if  del==-0 . 5 

[rmx, rmy, rmz] =rmhvector (rx, ry, rz,m-l) ; 

end 
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delta=rmag; 
delta2=delta/ 2  ; 

distl=sqrt ( (rmx-rnx) A2+ (rmy-rny) A2+ (rmz-rnz) A2) ; 
dist2=sqrt ( (rmx-rnx2) A2+ (rmy-rny2) A2+ (rmz-rnz2) A2) ; 
if  (distl<delta) & (dist2<delta) 

%  Self  term  calculation 


F=@ (s)  exp (-j *wk*sqrt ( (rmx-rnx-s*sxl) . A2+ (rmy-rny-s*syl ) . A2+ (rmz- 
rnz-s*szl) . A2+rad. A2) ) . / . . . 

sqrt ( (rmx-rnx-s*sxl) . A2+ (rmy-rny-s*syl ) . A2+ ( rmz-rnz - 
s*szl)  . A2+rad. A2) ; 

psi=quad (F, 0 . 0, delta) ; 


else 


%  Non-self  term  calculation  -  (Reduced  kernel  Approximation) 

F=@ (s)  exp (-j  *wk*sqrt ( (rmx-rnx-s*sxl)  . A2+ (rmy-rny-s*syl )  . A2+ (rmz- 
rnz-s*szl) . A2+rad. A2) ) . / . . . 

sqrt ( (rmx-rnx-s*sxl) . A2+ (rmy-rny-s*syl ) . A2+ ( rmz-rnz - 
s*szl)  . A2+rad. A2) ; 

psi=quad (F, 0 . 0, delta) ; 


end 


%  &&&&&&&&&  FUNCTION  "vecpot "  &&&&&&&&&&&&& 

function  psi=vecpot (rx, ry, rz, wk, rad, m, n, del) 

%psi  function  due  to  vector  potential 

rmx=rx (m) ; 
rmy=ry (m) ; 
rmz=rz (m) ; 
if  del==0 . 5 

rnx=rx (n) ; 
rny=ry (n) ; 
rnz=rz (n) ; 

[sxl,syl,szl, rmag] =sunit ( rx, ry, rz , n) ; 

end 

if  del==-0 . 5 

[rnx,  rny, rnz] =rmhvector (rx, ry, rz, n-1) ; 

[sxl,syl,szl, rmag] =sunit (rx, ry, rz, n-1) ; 

end 

delta2=rmag/ 2.0; 
if  m==n 
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F=@ ( s )  exp ( - j  *wk*sqrt ( (rmx-rnx-s*sxl )  .  A2+  ( rmy-rny-s*syl )  . A2+ ( rmz- 
rnz-s*szl) ,A2+rad.A2) )./... 

sqrt ( (rmx-rnx-s*sxl) . A2+ (rmy-rny-s*syl ) . A2+ (rmz-rnz- 
s*szl)  . A2+rad. A2) ; 

psi=quad (F, 0.0, delta2) ; 

%  Self  term  calculation 


else 

%  Non-self  term  calculation  (reduced  kernel) 

F=@ (s)  exp (-j  *wk*sqrt ( (rmx-rnx-s*sxl)  . A2+ (rmy-rny-s*syl )  . A2+ (rmz- 
rnz-s*szl) . A2+rad. A2) ) . / . . . 

sqrt ( (rmx-rnx-s*sxl) . A2+ (rmy-rny-s*syl ) . A2+ (rmz-rnz- 
s*szl) . A2+rad. A2) ; 

psi=quad (F, 0 . 0 , delta2 )  ; 

end 

%  &&&&&&&&&&&&&  FUNCTION  "efield"  &&&&&&&&&&&&&&&&&&&&&&& 
function  efield  = 

efieldnear (rx, ry, rz, xdd, ydd, zdd, solvector, wk, rad, factor,  delta,  wavelengt 
h, ix, iy, iz, nunkns) 

%  This  function  computes  the  near  electric  field  E  x  for  ix=l, iy=0, iz=0 

%  This  function  computes  the  near  electric  field  E  y  for  ix=0, iy=l, iz=0 

%  This  function  computes  the  near  electric  field  E  z  for  ix=0,  iy=0,  iz=l 

nunkns2=nunkns+2 ; 
mpl=nunkns2+l ; 

if  ix  ==  1 

%  ****  Compute  x-Component  of  Electric  Field  (E  x)  ********** 

rx (mpl ) =xdd-0 . 0005*wavelength; 
ry (mpl ) =ydd; 
rz (mpl ) =zdd; 

rx (mpl+1 ) =xdd; 
ry (mpl+1 ) =ydd; 
rz (mpl+1 ) =zdd; 

rx (mpl+2 ) =xdd+0 . 0005*wavelength; 
ry (mpl+2 ) =ydd; 
rz (mpl+2 ) =zdd; 

elseif  iy  ==  1 


34 


Distribution  A:  For  Public  Release 


k  k  k 


Compute  y-Component  of  Electric  Field 


(E_y) 


kkkkkkkkkk 


rx (mpl ) =xdd; 

ry (mpl ) =ydd-0 . 0005*wavelength; 
rz (mpl ) =zdd; 


rx (mpl+1) =xdd; 
ry (mpl+1 ) =ydd; 
rz (mpl+1 ) =zdd; 

rx (mpl+2 ) =xdd; 

ry (mpl+2 ) =ydd+0 . 0005*wavelength; 
rz (mpl+2 ) =zdd; 

else 

rx (mpl ) =xdd; 
ry (mpl ) =ydd; 

rz (mpl ) =zdd-0 . 0005*wave length; 


rx (mpl+1) =xdd; 
ry (mpl+1 ) =ydd; 
rz (mpl+1) =zdd; 

rx (mpl+2 ) =xdd; 
ry (mpl+2 ) =ydd; 

rz (mpl+2 ) =zdd+0 . 0005 *wave length; 

end 

%  ****Compute  z-Component  of  Electric  Field  (E  z)  ********** 


[rxl, ryl, rzl] =  rmhvector (rx, ry, rz, mpl+1) ; 

[ rx2 , ry2 , rz2 ] =  rmhvector ( rx, ry, rz , mpl ) ; 

dif fx=rxl-rx2 ; 
dif fy=ryl-ry2 ; 
dif f z=rzl-rz2  ; 

esum=0 . 0 ; 

for  n=l:nunkns 
npl=n+l ; 

%  Contribution  due  to  vector  potential 
psil=vecpot ( rx, ry, rz , wk, rad, mpl+1 , npl ,  0.5); 


35 


Distribution  A:  For  Public  Release 


psi2=vecpot ( rx, ry, rz , wk, rad, mpl  +  1 , npl  ,-0.5) ; 


%Contribution  due  to  scalar  potential 

psi3=scalarpot ( rx, ry, rz , wk, rad, mpl+1 ,+0.5, npl , npl+1 ) ; 
psi4=scalarpot ( rx, ry, rz , wk, rad, mpl+1 ,+0.5, npl-1 , npl ) ; 
psi5=scalarpot ( rx, ry, rz , wk, rad, mpl+1 ,-0.5, npl , npl+1 ) ; 
psi6=scalarpot (rx, ry, rz , wk, rad, mpl+1 ,-0.5, npl-1 , npl ) ; 

%  S  unit  vectors 

[sxl, syl, szl, rmag] =sunit ( rx, ry, rz , npl ) ; 

[sx2,  sy2,  sz2,  rmag] =sunit (rx, ry, rz, n) ; 

%  Dot  products 

dotl=psil* (dif fx*sxl+dif fy*syl+dif f z*szl ) ; 
dot2=psi2* (dif fx*sx2+dif fy*sy2+dif f z*sz2  )  ; 
dotprod=wkA2* (dotl+dot2)  ; 

mat rix=f actor* (dotprod-psi3/delta+psi4 /delta+psi5/delta- 
psi6/delta) ; 

esum=esum+matrix*solvector (n)  ; 
end  %  Loop  over  n;  Ex  Component 

efield  =  -esum/0 . 001 /sqrt (2 ) ;  %  x-comp  of  Electric  Field 
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